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Abstract 

The gravitational collapse of an infinite cylindrical thin shell of generic matter in an otherwise 
empty spacetime is considered. We show that geometries admitting two hypersurface orthogonal 
Killing vectors cannot contain trapped surfaces in the vacuum portion of spacetime causally avail- 
able to geodesic timelike observers. At asymptotic future null infinity, however, congruences of 
outgoing radial null geodesies become marginally trapped, due to convergence induced by shear 
caused by the interaction of a transverse wave component with the geodesies. The matter shell 
itself is shown to be always free of trapped surfaces, for this class of geometries. Finally, two 
simplified matter models are analytically examined. For one model, the weak energy condition is 
shown to be a necessary condition for collapse to halt; for the second case, it is a sufficient condition 
for collapse to be able to halt. 
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I. INTRODUCTION 



The issue of formation of horizons in generic gravitational collapse remains an outstand- 
ing problem of classical general relativity. The original concept of trapped surface, due to 
Penrose [|1|, defines it as a compact spatial two-surface S on which 6_6j^ > 0, where 6± 
are the expansions in the future-pointing null directions orthogonal to S. Trapped surfaces 
signal thus the boundary of a region where any initially expanding null congruence begins 
to converge; clearly, they define regions of 'no escape'. Assuming that cosmic censorship 
holds, the existence of trapped surfaces implies that of an event horizon, which contains the 
latter In spherical symmetry, both sufficient and necessary conditions for the occurrence 
of trapped surfaces have been obtained, which are essentially of the form m > r, where m 
and r refer to local (or quasi- local) definitions of mass and radius f^. 

Understandably, much less is known about horizon formation in non-spherical geometries. 
Arguably one of the strongest results to date is that of Shoen and Yau f^, who obtained 
a sufficient criterion for the formation of trapped surfaces in an arbitrary spacetime: for a 
given lower bound for the mass density, there is an upper bound for the matter radius, above 
which trapped surfaces will form. However, being a sufficient condition for the occurrence 
of trapped surfaces, this result cannot say anything about the conditions under which the 
collapsing spacetime fails to develop horizons, possibly leading to naked singularities. As 
with several other notable issues in relativity, a conjecture has been put forward regarding 
its solution: Thome's hoop conjecture states that "horizons form when and only when a mass 
M gets compacted into a region whose circumference in every direction is C < 47rM" As 
originally stated, the conjecture leaves ample room for different definitions of horizon, mass 
and circumference. In spite of this ambiguity, no known counter-example appears to exist: 
numerical simulations of prolate and oblate collapse 0, gravitational radiation emission 
in aspherical collapse 0, and analytical studies of prolate collapsing spheroids either 
confirmed or could not refute the conjecture. Detailed analyses of cylindrical pressureless 
dust 1^, 0, and counter-rotating dust collapse [ill] ^^^^ confirmed the non-occurrence of 
horizons. 

Cylindrical spacetimes (defined below) constitute an obvious class of spacetimes to 
study non-perturbative departures from spherical symmetry. Although such spacetimes 
do not model exactly the dynamical evolution — and inherent gravitational-wave emission — 
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of bounded bodies, they possess a non-trivial field content and constitute thus a valuable 
test-bed for numerical relativity |T^, quantum gravity |T3|, and for probing the cosmic cen- 
sorship and hoop conjectures [Q. Since several definitions of cylindrical symmetry exist, 
one must adopt one to work with. In this paper we shall consider cylindrical spacetimes 
defined by the existence of two commuting spacelike Killing vector fields (such that the 
orthogonal spacetime is integrable): one translational (dz) and the other with closed orbits 
((9^), where the azimuthal coordinate is to be identified at and 27r. In addition, we take 
dz and (9^ to be hypersurface orthogonal, which implies that the cylindrical waves in the 
vacuum regions admit only one polarization [|15], |16]. 

The purpose of this paper is twofold: (i) study the occurrence of trapped surfaces in 
dynamical cylindrical spacetimes with a thin shell of arbitrary matter, and (ii) motivate the 
notion that realistic (in a suitably defined sense) matter is necessary to prevent the formation 
of curvature singularities. Specifically, we consider an infinite cylindrical thin shell with the 
most general surface stress-energy tensor defined on it, in an otherwise vacuum spacetime. 
We examine three alternative criteria for trapped surfaces, and show that they can never 
form, regardless of the matter content for the shell. We also show that — for one class of 
matter models — the violation of the weak energy condition (WEC) implies that collapse 
cannot be halted by an outward pressure gradient, and therefore the formation of an infinite 
spindle-like singularity is inevitable. 

The paper is organized as follows: Section II derives the needed mathematical framework, 
and the vacuum and junction condition equations for the spacetime. Section III discusses 
trapped surfaces according to three alternative criteria, and shows that there are no trapped 
cylinders anywhere on the spacetime, except at the limiting case of asymptotic null infinity. 
This limiting behavior is then explained via the Newman- Penrose formalism. In Sec. IV, we 
consider two different classes of matter content, and show that, the WEC is either necessary 
for collapse to halt, or sufficient for it to have a chance of being halted. Section V concludes 
with a summary and discussion. 

Natural geometrized units, in which SnG = c = 1, are used throughout. 
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II. INFINITE CYLINDRICAL THIN SHELL OF GENERIC MATTER 



The complete four-dimensional spacetime consists of an interior vacuum region con- 
nected to an exterior vacuum by a three-dimensional thin shell S. The vacuum regions 
are characterized by the Einstein-Rosen metric ||T6| 



dsl = e2(^±-^±)(-rft| + drl) + e^'^^dz^ + rle-^^^^d^^, (2.1) 

where 7 = 7(t,r), ip = ip{t,r), and the coordinate systems {xj^} are adopted. On E there is 
a natural holonomic basis {e(a)} given by 

= (2-2) 

where {^", a = 0, 1, 2} are intrinsic coordinates on E. The induced three-metric on E is 
then 

7afc = ^^.^ef^^e^f,), (2.3) 

and it is the same on both sides of E, since the four- metric must be continuous across it. 
This leads to 

dsl = -dr^ + e'^^^^Uz^ + R''{T)e-^^^^^^d(l)^, (2.4) 



^ = Vi?2 + e2W5:-7±) = r/±. (2.5) 
dr 

The zz and (pep components of the metric junction condition [g^^] = g^^ — g'^^ = imply 
that the r coordinate is continuous across E. In the equations above, r is the proper time 
measured by an observer comoving with the shell, with four-velocity 

where R = dR/dr. The shell E is defined by 

$(a;^) =r-rE(t±) = r-i?(r) =0. (2.7) 
The spacelike unit normal to E is 

< = a-i(r)^f a±$ = R6^^ + r]±6i^, (2.8) 
where a = e'^^'^^~^^^r]^^ is a normalization factor. 



4 



The normal extrinsic curvature, i^^^, is |T7 



K. 



ab 



Its non-vanishing components are: 



d'^x" dx^ dx" 

+ r 

I in 



V± 



i? + i?(7±-^s) 



where D± = ti^V^ is the normal derivative with respect to S, and the object 



dt 



(2.9) 



(2.10) 
(2.11) 
(2.12) 



(2.13) 



was used. 

For the matter content on the shell, we consider the most general stress-energy tensor 
Sab defined on S, compatible with the metric ( ^.41) : 

(2.14) 



Sab = p6:6l + p.e^^H:5', + p^R^e-'^^Ht 



where p, p^, and p^ are the proper surface energy density, pressure along the z-direction, 
and surface stress, respectively, as measured by an observer comoving with the shell with 
four- velocity given by (p^). 



A. Vacuum field equations 



The vacuum field equations are (where the ± subscript has been dropped for simplicity) : 



7,r = r[itt)' + itrf]- 



(2.15) 
(2.16) 
(2.17) 



We are thus left with the axisymmetric wave equation ( p.l5|) , which is decoupled from Eqs. 
( p.l6| )-( p.l7| ). One first solves the wave equation with appropriate boundary data, and then 
solve for ( |2.16| )-( pTT7| ) by quadratures. Remarkably, these two equations are compatible 



because their integrabihty condition is precisely the wave equation ( |2.15| ). 
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The general (outgoing) solution to the wave equation ( |2.15|) is given by 

POO 

V'(t,r) = / A{uj)e-''^'H^^\ur)duj, 



(2.18) 



where A{uj) is a complex- valued function, and H, 



XI 



Jnix) + iYn(x) is a Hankel func- 



tion of the first kind, and Jn and Yn are Bessel functions of the first and second kind, 
respectively |jl8 . 



B. Junction conditions for the shell 



On the shell, Einstein's equations reduce to the Darmois-Israel junction conditions for 



the normal extrinsic curvature 17, IS 



[Ka,] = K^, - K 



ah 



-Sab + -S'^a 



The non-vanishing components of the above equation are zz, 
respectively: 

Iv] = -Rp, 



2 ' 



Rp 

V+V- 



{R-Rijj:), 



where 



(2.19) 

and TT, which yield, 

(2.20) 
(2.21) 
(2.22) 

(2.23) 



Using the field equations (|2.16|) - (|2.17|) , Eq. ( |2.22|) can be cast solely in terms of the wave-field 



(2.24) 



Equations (p. 201) and ( p.21| ) match the metric functions ip and 7 across S, and Eq. (p.24|) 
governs the motion of the shell. We note that, in the absence of a matter shell (p = 0), it 
immediately follows from Eqs. ( p.5|) and (|2.21|) that the 'time' coordinate t is continuous 
across S and hence it can be globally defined, as expected. 

An additional constraint is obtained by the conservation of energy- momentum on E: 



(3)v7 ob — _\f,t^ 7^ 



n 



the only non- vanishing component of which gives 

R 

P + {p + P4>)-5 + (Pz - P<t>)ips = [et)Ti,un 



R 



(2.25) 



(2.26) 
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III. TRAPPED SURFACES 



A. Proper circumference radius criterion 

We shall use the following operational definition of trapped surface, which is a natural 
extension of Penrose's for surfaces that are not necessarily compact nor two-dimensional: 
let S be any non-null surface and 6± the expansion in the future-oriented null directions 
normal to 5*, then 5* is said to be a trapped surface iff dj^6-\s > 0, where the inequality 
saturates for the marginally trapped case. In the present case, our surface S is that of an 
infinite cylinder with proper circumference radius TZ = re~^. 

Modulo the two-dimensional quotient space orthogonal to the symmetry axis, the 
Einstein-Rosen metric is conformal to Minkowski spacetime (which is why the four- 
dimensional wave equation is formally the same as that of fiat space) , and thus it suffices to 
introduce fiat spacetime null coordinates (where, for simplicity, the '±' subscript is omitted): 

u = t — r, V = t + r, (3.1) 

to study the expansion of 5* along normal null directions. We have then 

e± = {dt T dr)n. (3.2) 

Trapped surfaces form when 

e+e. = (7^,^)2 - {Urf > o. (3.3) 

1. Vacuum regions 

In the vacuum regions, the above condition reads 

r^Vj > (1 - rip^rY. (3.4) 

One can use the asymptotic properties of Bessel functions (cf. Appendix B) to show that, 
for ujr ^ 1, 

^ -^^^ ^ 3? / cuA(cu)e-*t"(*+")-3"/^l(ia;. (3.5) 
v'^ Jo 

Hence, for ur ^ 1, 

e+e_ ~ -1 + 2ri)^ti^^r < o. (3.6) 



For large radii, this regime covers 'almost' all the frequencies, whereas for small radii, the 
validity of the inequality above is restricted to high-frequency modes with uj ^ r^^ ^ 1. 

In general, it proves more convenient to look for outer marginally trapped surfaces 
(OMTS) — surfaces whose expansion vanishes along future-oriented null normal directions — 
since their existence is implied by that of trapped surfaces Hence, it suffices to show 
that no OMTS form, to prove that there are no trapped surfaces. The condition for OMTS 
is 

0+ = -i^,t - - + = 0. (3.7) 
r 

Partial differentiation of this equation with respect to t and r yields, respectively, 

i^,tr = i^,U, i^,tr = i^,rr + (3.8) 



Combining these two equations with the wave equation ( |2.15D leads to 



= l^,t = --■ (3.9) 

2 r 

Now, from Eq. ( |3.5| ), we know that these two equalities fail to hold for ur ^ 1. To show 
that the same happens for cur < 1, we note that 



/ x"^ H'^\x)dx < oo iff n > -1, 
Jo 



(3.10) 



from which it follows that, as ^ 0, A{uj) ~ uj"" with n > —1. Hence, for ur < 1, we have 

^.t-V'.r-^, (3.11) 

thereby showing that Eq. ( |3.9| ) also fails to hold for ur < 1. By construction, all of the 
above holds inside and outside the shell, and we conclude therefore that there are no trapped 
surfaces in neither vacuum region. 



2. Shell 



On the shell, the condition for OMTS is 

Since ri/R> 0, an obvious necessary condition for OMTS is 

|-V's>0. (3.13) 
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However, for a collapsing configuration the physical circumference radius must decrease with 
proper time: 

7^< 0^ < 0, (3.14) 

which contradicts the necessary condition for OMTS. Hence, a collapsing shell can never 
become trapped, irrespective of the details of the matter content. 



B. Specific area radius criterion 

An alternative definition of trapped cylinder has been recently proposed by Hayward |^ 



and relies on the specific area radius, which is constructed as follows. The two spacelike 
Killing vectors {dz, d^} have a well-defined invariant geometric meaning: 



n := \d^\ = (3.15) 

is the proper circumference radius, and 

I := m = ^/^~^'z (3.16) 

is the specific Killing length. One can thus define a specific area for the cylinder as 

A:=2M = 27rr, (3.17) 

where f = 7U is the specific area radius. A cylinder is said to be trapped, marginal, or 
untrapped if the vector A^^ := V^^f is timelike, null, or spacelike, respectively In our 
case, we have 

'jZ = re-^, l = e^, f = r. (3.18) 

From a geometrical viewpoint, this definition is perfectly acceptable, since, modulo invariant 
rotations, all of the spatial ('radial') submanifold orthogonal to the axis is covered by r G 
[0, +oo). However, from a physical standpoint, this definition is somewhat lacking, since the 
specific area radius is simply a coordinate radius, and as such changes under the rescaling 
2; — >■ as r — > a;~^/^r, unlike the proper circumference radius, which remains invariant 
under such rescaling. In addition, an external timelike observer can only measure proper 
circumferences, but not coordinate radii. Despite these shortcomings, and because of its 
useful geometrical meaning, we shall examine the formation of trapped surfaces according 
to it. 
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1. Vacuum regions 



In the vacuum regions, 

Nr = (3-19) 

and thus 

{N^N>^)^^^ = e2('^-^) > 0. (3.20) 

Hence there are no trapped surfaces in the vacuum regions provided the quantity ip — j does 
not diverge negatively (in which the case a marginally trapped surface would form). Since 
Qrr = ~gtt = e'^^"'~^\ it follows that for reg'u/ar spacetimes the inequality must be strict and 
thus there are no trapped surfaces. This is true for any cylindrical regular metric with two 
commuting hypersurface orthogonal Killing vectors. Accordingly, it holds inside and outside 
S. In the interior vacuum region, it is straightforward to show that regularity of the axis 
implies absence of trapped surfaces. Let x = |f^</>P = T^e~^'^; then the symmetry axis is 



regular iff the following local flatness condition is obeyed |]22 

linil^ = l, (3.21) 

which requires 

7(t,0) = 0, tp,t\r=o<^^ ^Ar=o<^- (3-22) 
Hence, |A^|vac > along the axis. 



2. Shell 



On the shell, 

iVj = Rr^'X + '^M' (3-23) 

which gives 

(N^m)^ = 7]-^e^^^-^^ ijf - R^) = 7]-^e^^^-^^ > 0, (3.24) 

where the junction condition (|2.5|) was used in the last equality. Since the metric must be 
regular on S, the inequality can never saturate, and thus S is free of trapped surfaces. 
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C. Geodesic null congruences 



A direct, physical way to probe the occurrence of trapped surfaces in the vacuum space- 
time(s), is to examine the divergence of a congruence of outgoing radial (i.e., orthogonal to 
the symmetry axis) null geodesies (ORNG). To examine the behavior along null directions, 
it is convenient to introduce retarded Bondi coordinates {u,r, z,(j)}, in which the Einstein- 
Rosen metric reads 

els'" = -e2(^-^) {du^ + 2dudr) + e^^dz^ + r^e-^^d(j)^. (3.25) 
In these coordinates, the vacuum field equations are 

2^^ur-^,r+r-\^,u-^,r) =0, (3.26) 

7,„ = 2rV',«(V^,« + ^,r), (3.27) 
l,r = r(^,.)'. (3.28) 

The null vector field tangent to the ORNG is given by 

k'' = ^ = X{u,r)6';, (3.29) 
dX 

where A is an affine parameter along the geodesies, and X > since the geodesies are 
outgoing (i.e., r is a monotonically increasing function of the affine parameter). The geodesic 
equation reads then 

rV.A;^ = X[X,r + 2X(7,, - tP^rM = 0, (3.30) 
which readily integrates to 

X{u,r)=f{u)e'^^-^\ (3.31) 

where f{u) > is an arbitrary function, fixed by the choice of X for a given value of the 
coordinate r. For such a geodesic congruence, the expansion scalar is 

e := V^A;^ = = (3.32) 

Clearly, as long as one restricts oneself to finite values of r, there are no trapped surfaces. 
There could be trapped surfaces at asymptotic future null infinity |2^ jT+if, for fixed u, 
lim^^+oo-^ ^ 0{r). From a physical viewpoint, the possible occurrence of marginally 
trapped surfaces at J'~^is not worrisome, in that geodesic timelike observers can never reach 
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[|2^. However, from an asymptotic infinity viewpoint, it is of interest to examine the 
limiting behavior of geodesic focusing in a radiative spacetime. The asymptotic behavior of 
9 at i7+is given by 



e 



hm 

r— >+co f 



u=uq 



liin ■^(^o) c2(V.(«o,r)-7(«o,r))_ 



(3.33) 



The behavior of the wave- field ip [which determines 7 via the field equations ( p.27| )-( p.28| )] 
at J7'"''was computed by Ashtekar, Bicak, and Schmidt [E3, 12^ , and yields 



ip{u, r) 



'r \/r 



^ 00 



u)r 



n=l 



l{u, r) = 70 - 2 / CqcIu - ^ 



bJu) 



n=l 



n + 1 ri+'' 



(3.34) 
(3.35) 



where Cq = dco/du, and the coefficients c„ and 6„ are determined from the initial Cauchy 
data. From Eqs. ( p.33|) -( p735l) it follows that 



ei 



0. 



(3.36) 



This limiting behavior reflects the fact that, near J'^the cylindrical waves behave locally 
like plane waves, which intersect the congruence of ORNG, thereby causing them to start 
(marginally) converging (cf. Fig. 1). 



Such behavior is best described using the Newman- Penrose formalism ^7j. The two 
relevant equations are (cf. Appendix C) 



De = e'^ + a(T + %o, 
Da = a{e + e) + ^^o), 



(3.37) 
(3.38) 



which govern the evolution of the expansion —6 and shear a of the ORNG, respectively. 
Near J'~^the congruence is orthogonally intersected by plane waves aligned with (cf. Fig. 



1). The existence of a transverse wave component in the n'^ direction |^ is signaled by a 
non- vanishing ^E'(o): 

*(o) = -C^pf^urmfl^m" ^ 0. (3.39) 

This non-zero \E'(o) implies — by Eq. ( 3.38| ) — that the congruence starts to shear, with the 
shear axis being determined by the polarization of the wave (i.e., the wave-fleld iJj). This 
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e=o 



Initially expanding 
geodesic congruence 



e>o 



Plane wavefront 



FIG. 1: Asymptotic behavior of an initially expanding future-oriented null geodesic congruence at 
. When the radius (proper or otherwise) of cylindrical waves tends to infinity, they become 
plane waves, whose wavefronts orthogonally intersect the null geodesic congruence. The transverse 
wave component in the direction induces shear on the congruence, which in turn makes it focus, 
by virtue of the Raychaudhuri-type equation (3.37). 



introduces a non-negative term 



^.2 



acy = ( V^,, - + — ) (3.40) 



into Eq. ( |3.37| ), thereby inducing the congruence to converge. 



D. C-energy 



The concept of C-energy was introduced by Thorne as a local definition of energy 
density per unit Killing-length z for cylindrical systems. Being a locally defined quantity in 
spacetimes with a translational Killing field, its relation to energy-like quantities defined in 
asymptotically flat spacetimes (e.g., the ADM or Bondi mass) becomes unclear In fact. 



as shown by Ashtekar and Varadarajan [pl^, 32 1, the Hamiltonian that generates asymptotic 



time translations at spatial infinity is not the C-energy, but a non-polynomial function of 
it, which is positive and bounded from above. Nevertheless, being a locally conserved and 
measurable quantity, C-energy remains a very useful tool for the analysis of cylindrical 
systems, and, unsurprisingly, it can be linked to the occurrence of conical singularities ||3^ 
and trapped surfaces pi] . 
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Following Thome we introduce the C-energy potential 



C(t,r):=i(l-e-^^), (3.41) 

which is proportional to the total C-energy contained inside a cylinder of specific area radius 
r per unit Killing-length z. The C-energy fiux vector is given by 

:= ^e'^^"^v^[ln(l-8C)],,^^ 

I 



2 

2 I 

p2(i/<-7) 

= (C,,5f-CA^), (3.42) 

r 

and obeys the local conservation law V^P^ = 0. An observer with 4- velocity measures a 
C-energy density E = P'^u^, and the C-energy fiux across a hypersurface S with spacelike 
normal (such that n'^u^ = 0) is JF = P^n^. 
From Eqs. (ICTp and ([CTD it follows that 

(iVMiV^)vac = e^^^ (1 - 8C) . (3.43) 
Hence, trapped surfaces form whenever 

^ > ^, (3.44) 

where the inequality saturates for the marginally trapped case. This is the cylindrical 
analogue of the well-known asymptotically fiat spherically symmetric condition E >l/2 for 



Tolman-Bondi spacetimes |Q, where E = mus/Rp is the Misner-Sharp mass per unit 



proper area radius. In the present case, C < 1/8 and the vacuum regions are untrapped. 

IV. SINGULARITIES 

The issue of singularity formation in cylindrical collapse of a thin shell has been addressed 
for the particular cases of pressureless dust and counter- rotating dust (i.e., with a non- 
vanishing tangential pressure component). The pure dust case was analyzed by Thorne 
and Echeverria []TU[, who showed that complete collapse is inevitable, and an infinite spindle- 



like singularity must thus form as a result. The case of counter-rotating dust coupled to 



Einstein- Rosen waves was addressed by Apostolatos and Thorne |jTT|, who resorted to se 



quences of momentarily static radiation-free (MSRF) configurations and C-energy balance 
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arguments to conclude that an arbitrarily small amount of counter-rotation (and, very likely, 
also of net rotation) suffices to halt collapse, thereby precluding the formation of spindle 
singularities. The counter- rotating case was recently revisited by Pereira and Wang p6| . 



who considered a flat interior and an outgoing null fluid exterior (as opposed to a 'realistic' 
vacuum containing cylindrical waves); by means of further assumptions for the matter con- 
tent, they obtained a simple solvable model which admits initial data leading to a spindle 
singularity, thus in apparent contradiction with the earlier results. However, the model is 
highly contrived, and the assumptions made render it unphysical, since the dynamics of the 
exterior null fluid can only be obtained a posteriori, after a particular solution for the motion 
of the shell is derived. It is therefore conceivable that a more realistic model, wherein the 
emission of gravitational radiation arises dynamically from the inward accelerated motion of 



the shell, will qualitatively agree with the results of |11 



We shall conjecture here that realistic matter is required to prevent the occurrence of 
singularities. To support this conjecture, we examine below two classes of matter models 
for which the WEC is either necessary for collapse to halt, or sufficient for it to be able to 
halt. 

A. A = 

We consider here the case A = p = — pz- This case includes the Pereira- Wang 
matter model {pz = 0, p = p^), and a subcase of the Apostolatos-Thorne model {pz = 0, 
but not necessarily p = p^). Then, the junction conditions ( |2.20| ) and ( ^.241 ) simplify to 

= 0, (4.1) 

R = R^^ + V+V-^ - Rii^l + {D±iJ-f}, (4.2) 

and the local conservation equation on the shell becomes 

p V p J R 

This equation can be easily integrated by assuming a functional relation between p and p^, 
which, for simplicity, we take to be of barotropic form: p^ = ap, where the real constant 
< 1, to preserve causality. This gives then 

P = ^^' (4-4) 
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where A > is an integration constant, which is dimensionful and scales as [A] ~ i?" ~ L", 
such that [p] ~ L^^. This constant A turns out to be constrained by the fact that there is an 
upper hmit for the total rest mass per unit proper length, which arises from the requirement 
of radial nonclosure of space. Specifically, the space around a given cylindrical configuration 
is radially non-closed iff 

2nRaAT < \, (4.5) 

where ctat = p/ (8vr) is the total rest mass per unit proper area (the factor of Svr arises 
because we work with units in which SttG = 1, whereas Apostolatos and Thorne use G = 1). 
In our case, the total rest mass per unit proper length is the dimensionless quantity 

rrip = ^Rp = ^Ae'^^i?-", (4.6) 

from which it follows that we must have 

Rp < 1, (4.7) 



or A < e ^^R°', if ( [4.4|) holds. (We note that this constraint on the rest mass per unit proper 



length plays no role in the subsequent assertion that the WEC is a necessary condition for 
singularities not to form.) 

We now look for local extrema of R, and show that the obedience of the WEC is necessary 
for collapse to be halted. Local extrema of R occur at r = Te, such that Re = and Re 7^ 0. 
From Eq. ([4. 21) we have then 



Re = -V-iP + Pz) [^-^pj' ^i^s + iD±^~y}, (4.8) 

where the junction condition ( p.21|) was used, and the right-hand-side is to be evaluated at 
T = Te- Since the last summand is manifestly negative and //_ > 0, a necessary condition 
for collapse to be halted (i.e., for Re > 0) is 



Rp 

By Eqs. (|2.5| ) and (|2.21|) , the second term is always positive for p > 0. It then follows that 



we must have p + > if collapse is to be halted. This means that the WEC (which 
requires p > and p + Pi > Vz) must be satisfied in order for R to admit a real positive 
minimum; in addition, since p = Pif, — Pz hypothesis, we must also have p^f, > 0. If 
Pz < —p and/or p^ < 0, then condition ( [4.9|) is violated and there are no local minima of R, 
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with collapse inevitably proceeding to i? = 0, wherein a spindle-like singularity forms. We 
remark that this case requires at least one of the principal pressures to be non-vanishing, 
since = = ^ p = 0. 



B. A = /5 

This case corresponds to isotropic pressure, Is. = p <^ pz = P4>- The acceleration equation 
at the local extrema reads 



Re = V-[ V+f-p + ^±^-J - ^ - mi + {D^^_Y}. (4.10) 



At the extrema, = 0, and, since A = p, it follows from the junction conditions ( p. 201) . 
( ^ ) that 



Equation (^.lOj) can then be rewritten as 



The last summand is manifestly negative, so a necessary condition for > is 

vA + \v->0, (4.13) 
p 4 

where all quantities are evaluated at r = Tc. Clearly, if the WEC holds, the condition above 
is satisfied. Condition (|4.13| ) defines a negative lower bound for p^, 



P, > -\p'-, (4.14) 

which would be incompatible with a violation of the WEC (i.e., with < — p) iff 4?7+ > rj^. 
By the junction condition ( |2.21|) , this is equivalent to r]^ > {Rp)/3, where rj'^ = e^^~"'+ . Now, 
for a given r^. = -R(r^,), we can rescale the Killing coordinate z, such that = ipT:{R*) = 0; 
doing this for Rg, the inequality t]'^ > {Rp)/3 becomes 

7^ < -ln(i?p/3). (4.15) 

To examine whether this can at all be verified, we will assume that, in addition to being 
momentarily static (which, by definition, is the case at the extrema), the configuration is 
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also radiation- free at r = Xg, i.e., that ip^t = 4',tt = at t±{Te). With these ansatze, the 
wave-equation (p.l5| ) has the trivial solution: 



^± = V's-A;±ln(r/i?), (4.16) 
where k± is a constant. Then, Eqs. ( ^.16[ ),( P?T7D yield 

7± = 7± + 4M^/^)- (4.17) 



The regularity conditions ( |3.22 ) imply /c_ = 7^ = 0, thereby fixing the interior solution as 



ip- = ip-E and 7_ = 0. From Eqs. (|2.5| ) and ( |2.21|) , together with the solution for 7_, we get 



7; = -ln(l -i?p). (4.18) 

Hence, condition ( |4.15|) would be automatically satisfied (i.e., the WEC would be a necessary 
condition for collapse to halt) iff 

Rp<l- (4.19) 

This inequality is stronger than the one imposed by radial nonclosure of space [Rp < 1), and 
as such may or may not hold at r = Te. For this particular case, the approach adopted here 
cannot show that a violation of the WEC contradicts condition ( [4.13| ) without additional 
assumptions for the initial data (this does not, of course, imply that a violation of the WEC 
necessarily allows for collapse to be halted). This case is nevertheless instructive in that 
it (i) shows that the WEC is a sufficient condition for collapse to have a chance of being 
halted, and (ii) provides a self-consistency check: if the WEC holds, then the requirement 
of radial nonclosure of space is automatically satisfied. 



V. CONCLUSIONS 



We studied in detail the formation of trapped surfaces in polarized cylindrical spacetimes 
with a thin shell of arbitrary matter. We showed that, for arbitrarily large r, no trapped 
surfaces form, either in the vacuum regions, or on the shell, regardless of the matter content 
of the latter. We studied the limiting behavior of the expansion of a congruence of future- 
directed outgoing radial null geodesies, and found that they become marginally trapped 
exactly at future null infinity. This limiting behavior is due to the interaction between 
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the null geodesies and an orthogonal (locally) plane wave, whose non-vanishing transverse 
component induces the congruence to start focusing, via a Raychaudhuri-type equation. 

The fact that no trapped surfaces form irrespective of the matter content for the shell, 
could suggest that this is perhaps an artifact of the simple geometry adopted, in which case 
it would not necessarily hold in a fully generic cylindrical spacetime, with non-hypersurface- 
orthogonal Killing vectors (and whose vacuum regions, because of this, admit two polariza- 



tions). However, preliminary calculations |^9| indicate that this may not be the case, thus 
in accord with the spirit of the hoop conjecture. 

Even for the simple case of irrotacional dust {pz = p,j> = 0) , an analytical solution for the 
shell motion in closed form appears impossible to obtain, with the only available fully dy- 
namical solution being numeric Since there is no pressure, total implosion is inevitable, 
and one can thus derive analytical approximations for the late stages of collapse, when the 
dynamics is highly relativistic, and the shell asymptotes null dust. Such approximations are 
valuable for an analytical treatment of the singularity and its structure. In the presence of 
pressures, however, total implosion is not guaranteed (in fact it was shown not to occur, in 
the case of couter-rotating dust |Tl|) and, accordingly, one cannot use any approximation for 
the late stages, since one cannot know a priori whether and/or when collapse is halted. To 
be able to make statements about the endstate of collapse without numerical aid one must 
thus resort to specific classes of models and look for simplifying (but admissable) ansatze. 

To examine the effects of matter on the occurrence of singularities, we considered here 
two classes of matter models, which generalize and complement existing examples in the 
literature. By using the junction conditions on the shell, together with the assumption of 
positive energy density, we showed that, for the first class the WEC is a necessary condition 
for collapse to halt, and a sufficient condition for collapse to be able to halt in the second 
matter model. Based on these results, and those of pressureless and counter-rotating dust, we 
have put forward the mild conjecture that realistic matter is required to prevent singularity 
formation. In the conjecture, realistic is meant to imply that the WEC holds and at least 
one of the principal pressures is non-vanishing. 

The inclusion of another polarization makes the model mathematically more complicated, 
but the physics appears to remain largely unchanged (in particular, no singularities seem to 
form) . Accordingly, we expect the trapped surface results presented here to qualitatively 
extend to the unpolarized well as the conditions regarding the shell dynamics. Work 
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in this direction is currently underway 
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APPENDIX A: NON-EXISTENCE OF A GLOBAL FOUR-DIMENSIONAL CO- 
ORDINATE SYSTEM 

Let us consider the stress-energy tensor associated with S expressed as a four- dimensional 
distribution: 

T^" := S'^'e'l^fl.^lam). (Al) 

Its components are: 

The independent non-vanishing components of G^y are: 

Gtt = Grr = - + ^"S,r, (A3) 

Gtr = -2V^,,7/;,, + r-S,i, (A4) 

G.. = e-27+^'^ [2{^,u - ^,rr - r- V.r-) + 1,rr ' 1 ,tt ' + ] , (A5) 

G<^0 = e-^^r^ (7,.. - 7,« - V^J + ^J) • (A6) 

Einstein's equations give then 

-i^l - ^2 + = Tu = Trr, (AT) 

7,t = 2rV',tV,r, (A8) 

1p,tt - ■ip,rr - r"V,r = ^{Pz - P^)S{<^). (A9) 

21] 



In the vacuum regions these equations reduce to the set (p.l5|) -( ^rT7D , as expected. On 



the shell, however, they imply = =^ rf = E?, which is only true if [cf. Eq. ( p.5|) ] 
7^ — i> -|-oo, but this in turn renders the four-metric singular on S. This simply means 
that one cannot use a four-dimensional formulation for the whole spacetime, since the (four- 
dimensional) Einstein equation = T^r only holds in the vacuum regions. One must 
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therefore adopt two different coordinate systems, for the inner and outer vacuum regions, 
and match them across S via the Darmois-Israel junction conditions. 



APPENDIX B: ELEMENTARY PROPERTIES OF BESSEL FUNCTIONS 

Consider the integral 

/•oo 

J(t,r)=/ A{uj)e-''^'H^^\ujr)duj, (Bl) 
Jo 

where A{uj) is a complex-valued function, and Hn\x) := Jn{x) +iYn{x) is a Hankel function 
of the first kind, and J„ and are Bessel functions of the first and second kind, respectively. 
Partial differentiation of I yields 

POO 

It = - iuoA{uo)e'"^^H^^\ujr)duo, (B2) 
Jo 

POO 

I,r = - ujA{uj)e-''^'H[^\Lur)duj. (B3) 
Jo 



Now, for z ^ 1, we have 



i7(i)(^)^./Ae^[-f("+f)], (B4) 

V TIZ 



from which it follows that, for tur ^ 1, 



^ r .f^A{uy^-^^-'^-^Mu, (B5) 
Jo V Txr 

Jo V vrr 

^ , r ./5A(cu)e^[-(^-*)-f]f/a;. (B6) 
Jo V Trr 



Hence, I^r ~ — for cur ^ 1, thereby justifying Eq. (|3] 

APPENDIX C: NEWMAN-PENROSE FORMALISM FOR THE ANALYSIS OF 
LIMITING GEODESIC BEHAVIOR AT J+ 



Here we present a summary of the Newman- Penrose (NP) formalism [27|, which is par 



ticularly useful for a geometrical analysis of the behavior of null congruences. We introduce 
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a null tetrad, given by two real null vectors P and n'^, and a complex conjugate pair m'^ and 
m^, defined such that (i) their only non- vanishing inner products are 

l^n" = -m^m'^ = 1, (CI) 

and (ii) the metric completness relation holds: 

= lixriv + riijl^ - m^rhy - fh^rriy. (C2) 

With these definitions, the Ricci and Weyl tensors can be naturally decomposed in terms 
of their tetrad components. The Ricci tensor can be decomposed into a scalar component 
(which gives the curvature scalar) and a Hermitian (3 x 3) matrix <^ab which represents 
the trace- free part of the Ricci tensor and satisfies ^ab — ^ba- The ten independent 
components of the Weyl tensor Cap^i,, which represent the ten degrees of freedom of the 
gravitational field, are conveniently expressed as five complex NP scalars. The one relevant 
for our purposes is 

*(o) = -C^p^XmH^^m^ (C3) 

Each NP scalar has a distinct physical interpretation, associated with the presence of gravi- 
tational waves along the null directions and n^] ^'(o) denotes a transverse wave component 
in the direction. 

The behavior of geodesic null congruences (among other things) can be conveniently 
described in terms of NP spin coefficients, which are complex linear combinations of the 
Ricci rotation coefficients associated with the null tetrad. In the present case, the relevant 
spin coefficients are 

e = VJ^.m^fff, (C4) 
(7 = VJf,mV, (C5) 
k = VJi^m^r (C6) 

They measure the expansion, shear, and deviations from (null) geodesic motion for null rays 
along l^, respectively. The evolution of these quantities is given by directional derivatives 
in the directions of the four tetrad vectors: 

D = FV^, A = n^V^, 6 = m'^V^, 6 = m^V^, (C7) 
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from which one obtains 



D9-6k = e^ + aa-kr- k{3a + p - n) + $00, (C8) 
Da -5k = a{9 + 9) + a{e - e) - k{7r - n + a + 3/?) + ^(0). (C9) 

These equations govern the evolution of the expansion and shear along the null direction. 

The natural null tetrad for the Einstein- Rosen metric ( |2.1| ) is (where the metric signature 
has been changed to —2, to conform with the original NP construction): 

I, = e'^^-^\5l + 5;), (CIO) 
n, = liSl-6;), (Cll) 

m, = --Le^(z<5^ + re-2^5j). (C12) 
The non-vanishing Ricci components are: 

$00, '^02, '^'11, $22- (C13) 

The non-vanishing NP scalars are: 

^^(0), ^(2), *(4). (CM) 

Finally, the non- vanishing NP spin coefficients are (modulo complex conjugation): a, A, 6, n, 
and e. The relevant ones are: 

a = a = ip^t -■ip,r + (C15) 

2r 

e = e = -l, (CI6) 

e = e = i!.r-i)^t- 7,r + 7,i- (C17) 
The relevant evolution equations become then 

DO = e^ + aa + $00, (C18) 
Da = a{e + e) + ^(^oy (C19) 

All of the NP objects were computed with the aid of the algebraic package GRTensorll 
for Maple V. Most of the non-zero objects are 'too large' to be explicitly written in any 
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convenient way. Since only the fact that they are non-vanishing matters for our purposes 
(specifically, ^'(o) 7^ 0), we omit here their explicit form. 
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